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The optical properties of graphene are made unique by the linear band struc- 
ture and the vanishing density of states at the Dirac point [TT-llOj . It has been 
proposed that even in the absence of a semiconducting bandgap, a relaxation 
bottleneck at the Dirac point may allow for population inversion and lasing at 
arbitrarily long wavelengths |4H6]. Furthermore, efficient carrier multiplication 
by impact ionization has been discussed in the context of light harvesting ap- 
plications J7], [8]. However, all these effects are difficult to test quantitatively 
by measuring the transient optical properties alone, as these only indirectly re- 
flect the energy and momentum dependent carrier distributions. Here, we use 
time- and angle-resolved photoemission spectroscopy with femtosecond extreme 
ultra-violet (EUV) pulses at 31.5 eV photon energy to directly probe the non- 
equilibrium response of Dirac electrons near the K-point of the Brillouin zone. 
In lightly hole-doped epitaxial graphene samples | |11|, 112] . we explore excita- 
tion in the mid- and near-infrared, both below and above the minimum photon 
energy for direct interband transitions. While excitation in the mid-infrared 
results only in heating of the equilibrium carrier distribution, interband exci- 
tations give rise to population inversion, suggesting that terahertz lasing may 
be possible. However, in neither excitation regime do we find indication for 
carrier multiplication, questioning the applicability of graphene for light har- 
vesting. Time-resolved photoemission spectroscopy in the EUV emerges as the 
technique of choice to assess the suitability of new materials for optoelectronics, 
providing quantitatively accurate measurements of non-equilibrium carriers at 
all energies and wave vectors. 

The non-equilibrium dynamics of Dirac electrons in graphene have so far been primarily 
investigated with time-resolved optical techniques [5l [T3l - H8] . yielding precise estimates of the 
relaxation rates after photoexcitation. However, as optical probing only indirectly measures 
energy- and momentum-dependent carrier distributions, controversies have arisen over the 
underlying physics. Crucially, it has not yet been established if the distribution of hot 
carriers at early times should be described by a single [16], or by two distinct Fermi-Dirac 
(FD) distributions with different chemical potentials for valence and conduction band [5j 
TBI Ull US]- Clarification of this issue is essential for assessing the potential of graphene 
for optoelectronics as it defines whether charge carriers in graphene exhibit a metallic or 
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FIG. 1: (a) Honeycomb lattice with hexagonal Brillouin zone. The two triangular sublattices 
are shown in red and blue, respectively, (b) Equilibrium band structure measured with linearly- 
polarized synchrotron radiation {fojj = 30 eV, left) and He II radiation (fouj = 41 eV, right). The 
color scale is linear with violet (red) corresponding to high (low) photocurrent. The inset shows 
the cut through the Brillouin zone along which the photoemission data has been measured. Along 
the TK-direction the photocurrent is suppressed for one of the two 7r-bands |22j . The graphene 
layer is slightly hole-doped so that the chemical potential \i e is ~200meV below the Dirac point. 
Dashed black lines represent the dispersion obtained from a tight-binding model |20j. 



a semiconducting response to photoexcitation. Furthermore, direct experimental proof of 
the existence of carrier multiplication [TJ |8] is lacking. Here, we resolve these questions by 
directly mapping the transient occupation of electronic states as a function of momentum 
and energy with photoemission spectroscopy for different time delays after photoexcitation. 

The honeycomb lattice of graphene, together with a sketch of the two-dimensional Bril- 
louin zone is shown in Fig. [T^i. We studied a quasi-freestanding epitaxial graphene monolayer 
grown on SiC(OOOl) [TT| [1~2] . which was first characterized by static angle-resolved photoe- 
mission spectroscopy (ARPES). Figure [lb shows photoelectron energy and momentum maps 
along and perpendicular to the TK-direction, evidencing the well-known conical band struc- 
ture [201 EH- For measurements perpendicular to the TK-direction, extrapolation of the 
bands based on a tight-binding model (dashed lines) [20] shows that the graphene layer is 
hole-doped, with the Dirac point ~200meV above the static chemical potential \x e . For 
measurements along TK, only one of the two 7r-bands is visible due to photoemission matrix 
element effects [22J. 

The non-equilibrium response of this sample was studied for excitation below and above 
the minimum photon energy for interband excitations (2/i e = 400 me V). In this way, both 
metallic (ftcjpump = 300 me V) and semiconducting (foupump = 950 me V) optical properties 
were explored. To probe the dynamics of the 7r-bands at the K-point of the hexagonal Bril- 
louin zone (fen = 1.7A -1 , requiring hw > 16 eV), 30 femtosecond (fs) EUV pulses from a High 
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FIG. 2: Free carrier absorption regime for excitation at hco 



l>U mp = 300 meV: (a) Cartoon of the 
excitation mechanism. For hui pump < 2|// e | direct interband transitions are impossible. The Fermi 
liquid sloshes around the Dirac cone following the electric field E(t) of the laser pulse (light grey 
line) like water in a shaking glass. Electron-electron and electron-phonon scattering events sub- 
sequently establish a thermal Fermi-Dirac distribution. Within the time resolution of the present 
investigation the coherent dynamics at earliest times could not be resolved, (b) Snapshots of the 
band structure (smoothed) close to the K-point for different pump-probe delays: negative time 
delays (left), close to zero delay (+30 fs, middle), t = lps (right). The sample is excited at 
fofpump = 300 meV with a nuence of F = 0.8mJ/cm 2 . Photoelectrons are ejected using photons 
of energy ^Wp ro be = 31.5 eV. The measurements were done along the TK-direction (inset), where 
the photocurrent is suppressed for one of the two 7r-bands [22]. Dashed black lines represent the 
dispersion obtained from a tight-binding model 



Harmonic Generation (HHG) source were used J23H25]- A time preserving monochromator 
selected a single harmonic at hu pvo b c = 31.5 eV with energy resolution of 130 meV [26J. 

Figure [2^i shows a cartoon of the expected excitation mechanism for photon energies 
^pump < 2|/x e |, for which only free carrier absorption is possible. For an electromagnetic 
pulse with vector potential A(t) = J Edt, the Dirac electrons are accelerated and decelerated 
according to the interaction Hamiltonian H = ^*(p ~ qA) 2 . Scattering events with other 
electrons or with the lattice cause rapid thermalization into a hot carrier distribution. 

A set of experimental snapshots of the Dirac cone is plotted in Fig. ^p for different 
pump-probe delays after excitation at hu pump = 300 meV with a nuence of F — 0.8 mJ/cm 2 . 
Within the time resolution of our measurement we did not resolve non-thermal energy distri- 
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FIG. 3: Metallic response for excitation at hujp Urap = 300 meV and F = 0.8mJ/cm 2 : (a) Energy- 
distribution curves (EDCs) integrated over the momentum range displayed in Fig. ^p as a function 
of delay in a two-dimensional color plot, (b) Selected EDCs extracted along the dashed lines in 
panel (a) together with Fermi-Dirac distribution fits that are shown as solid black lines. Chemical 
potential (c) and electronic temperature (d) as a function of delay. The solid black line in (d) is 
a fit including an error-function for the rise time and three exponential decays: T\ = (13 ± 4) fs, 
7"2 = (80 ± 10) fs and T3 = (0.86 ± 0.05) ps, consistent with the results of all-optical measurements, 
(e) The fluence dependence of the maximum electronic temperature follows the F 1 ' 3 behavior 
(black line) expected for free carrier absorption. 



butions for the Dirac electrons, indicating very efficient thermalization. Rather, we observed 
instantaneous broadening of the carrier distribution through the Dirac point (Fig. [2b middle 
panel) and subsequent relaxation within ~1 picosecond (ps, Fig. [2b right panel). 

A quantitative analysis of these dynamics is shown in Fig. pj] In panel (a) we display 
energy distribution curves (EDCs) that have been integrated at each time delay over the 
momentum range depicted in Fig. [2]d. Zero delay corresponds to the maximum pump-probe 
signal. The energy dependence of the EDCs is well fitted at all times by the FD distribution 
/(.£/, fj, e , T e ) = ( exp f fc ~y e ) + 1 ) ? where E is the energy, /i e is the chemical potential, k B 
is the Boltzmann constant, and T e is the electronic temperature. The broadening of the 
Fermi edge immediately after absorption of the pump pulse reflects an increase in electronic 
temperature T e that subsequently cools down to room temperature within ~lps. Figure 
[3b displays EDCs for selected time delays [-0.3 ps (blue), ps (red) and +0.5 ps (green)] 



together with FD fits (black lines). Fit results for the chemical potential \x e and the electronic 
temperature T e are shown in Figs. [3b and d, respectively. 

The chemical potential initially increases by approximately 40meV, recovering as the 
carriers cool. The initial increase and fast decay of the chemical potential is well understood 
as a consequence of the changes in electronic temperature (Fig. 3m, as A/i = ^(/c^Te) 2 /^ 
[27] . where Ep is the Fermi energy. 

The time-dependent electronic temperature T e (t) is shown in Fig. |3]i. T e is observed 
to instantly rise to nearly 3000 K, before decaying with a triple exponential law (black 
line in Fig. [3J1), in excellent agreement with previous all-optical measurements. These 
measurements have attributed the decay channels to electron-electron scattering (n) [T31IT4"] . 
emission of optical phonons (r 2 ) [ToT - TITl 128] . and further cooling through the decay of optical 
into acoustic phonons (T3) [181 1291432] . More recently, the third relaxation time t% has been 
ascribed to the direct coupling between electrons and acoustic phonons in the presence of 
lattice defects [331 El]- Our data is consistent with both interpretations of r 3 . A validation 
of the decay law would require temperature-dependent measurements that are beyond the 
scope of his paper. 

The fluence dependence of the peak electronic temperature, which further substantiates 
the physical picture discussed above, is shown in Fig. [3^. The increase in kinetic energy 
of the Dirac carriers can be obtained from AE(T e ) = J DOS(E)f(E,fj, e ,T e )EdE, where 
DOS(E) oc \E\ is the linear density of states of graphene. For simplicity, we choose fJ, e (T e ) = 
and obtain AE ex T e 3 , where AE is equivalent to the absorbed fluence F in mJ/cm 2 . The 
experimental fluence dependence of T e , displayed in Fig. |3p, follows the expected F 1//3 
behavior (black line). 

We next report results for direct interband excitation at hoj pump = 950 meV in Figs. [4] and 
[5] Figure [4^i displays a schematic of the expected interband excitation mechanism, which, in 
contrast to the free carrier absorption discussed above, is driven by the intensity envelope I(t) 
of the laser pulse (red line) rather than the electric field. During interband excitation, a non- 
thermal distribution of excited electron- hole pairs is generated at the earliest time delays, 
relaxing into a quasi-equilibrium state with a temperature higher than that of the lattice. 
Importantly, depending on the relative strength of interband versus intraband scattering, a 
single metallic distribution like the one of Fig. [3]d, or two distinct distributions for valence 
and conduction band are attained. Only the latter scenario would make light amplification 
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FIG. 4: Direct interband transition regime for excitation at hu pnmp = 950meV: (a) Cartoon 
of the excitation mechanism. The red line represents the intensity envelope of the pulse. If 
liWp Um p > |2/i e |, the graphene layer is excited via direct interband transitions that leave the sys- 
tem with a highly non-thermal electron distribution that cannot be observed within the temporal 
resolution of the present experiment. Within ~ 10 femtoseconds the electrons thermalize and form 
separate Fermi-Dirac distributions for the valence and conduction band. After some hundreds 
of femtoseconds the separate chemical potentials merge, resulting in a single Fermi-Dirac distri- 
bution for valence and conduction band, (b) Snapshots of the band structure (smoothed) for 
excitation at huj purap = 950 meV with a fluence of F = 4.6mJ/cm 2 and a probe photon energy 
of ^Wprobe = 31.5 eV. The measurements were done along the TK-direction (inset). Dashed black 
lines represent the dispersion obtained from a tight-binding model 



possible. 

In Fig. [4b we present snapshots of the Dirac carrier distributions for different time delays 



after excitation at hu 



pump 



950 meV and a fluence of F = 4.6mJ/cm 2 . The excitation 
fluence was chosen according to experimental estimates of the threshold for optical gain in 
Ref. 0. 

Figures [5^ and [5b report momentum- integrated EDCs as a function of pump-probe delay. 
In contrast to Fig. [3b, the carrier distributions measured immediately after excitation are 
best fitted by two separate FD distributions for valence and conduction band, indicating 
population inversion [5]. The fit results for the chemical potential /x e and the electronic 
temperature T e are displayed in Figs. [5b and [5]i, respectively. The two distributions merge 
within ~130fs, and a single FD distribution is attained, promoted by scattering of charge 
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FIG. 5: Population inversion for excitation at foj pump = 950 meV and F = 4.6mJ/cm 2 : (a) 
Momentum-integrated energy distribution curves (EDCs) as a function of delay, (b) EDCs for 
three selected delays extracted along the dashed lines in (a) together with Fermi-Dirac distribution 
fits. The EDCs close to time zero are best fitted with two distinct Fermi-Dirac distributions for 
electrons and holes (dashed black lines). Chemical potential (c) and electronic temperature (d) 
as a function of delay. The red-shaded area indicates the time interval over which two distinct 
Fermi-Dirac distributions persist. The solid black line in panel (d) represents a fit of the black 
data points with a fixed rise time of 50 fs and a double exponential decay. The decay times are 
indicated in the figure. 



carriers across the Dirac point. As amplification is only possible for the given time window, 
our measurement sets a quantitative boundary to be used as a benchmark for any laser 
design. 

Our measurements also provide a direct assessment of the possibility of carrier multipli- 
cation in graphene [7J HJ. Figure |6] displays the time-dependent evolution of the pump-probe 
signal for the two cases of free carrier absorption (right panel) and interband excitation 
(left panel). Specifically, Fig. [6^, shows the pump-induced changes of the photocurrent 
(directly proportional to the number of electrons iV at energy E) as a function of energy 
and time delay. Blue (red) corresponds to a gain (loss) in photocurrent with respect to 
negative time delays. In Fig. [6)3, the time derivative of the pump-probe signal of panel 
(a) is displayed, directly revealing the electron-hole pair generation/recombination rates. 
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FIG. 6: No indication of carrier multiplication for direct interband transitions (left column) and 
free carrier absorption (right column): The pump-probe signal in panel (a) is given by the difference 
between the photocurrent shown in Figs. [3^, and [5^,, respectively, and the photocurrent at negative 
time delays. Panel (b) shows the time derivative of the pump-probe signal, which is equivalent to 
the electron- hole pair generation/recombination rate. Electron-hole pairs are generated only for a 
short time interval in the presence of the pump pulse. The data has been smoothed. 

The generation/recombination rate dN(E > 0,t)/dt is positive only for a few tens of fem- 
toseconds. We conclude that electron-hole pairs are generated only while the pump pulse is 
present, whereas no spontaneous carrier multiplication is found. 

The absence of carrier multiplication is well understood when analyzing the fits of Fig. 
[5j Indeed, in the presence of two distinct FD distributions, impact ionization [7J |8] can only 
create additional electron-hole pairs if the temperature of the Dirac carriers ksT e exceeds the 
sum of the electron and hole chemical potentials |/i e | + |/^| (Fig. [7k). This would require 
temperatures above ~8000K, significantly higher than the peak temperature of ~4000K 
estimated in Fig. [5] Instead, the inverse scattering process (Auger heating) dominates, 
rapidly merging the separate chemical potentials. For a single FD distribution (Fig. Up) 
attained ~130fs after direct interband transitions or at all times for free carrier absorption, 
the probabilities for impact ionization and Auger heating are similar, again preventing carrier 
multiplication. The absence of carrier multiplication and the short lifetime of electron-hole 
pairs discussed above raises serious doubts about the suitability of graphene for efficient 
light harvesting for the present excitation regime. One may speculate that more favorable 
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FIG. 7: Competing relaxation channels in graphene: For a population- inverted state with separate 
chemical potentials for valence and conduction band (a) impact ionization (green arrows) can only 
create additional electron-hole pairs if ksT e ~ |// e | + |/z/i|, which requires temperatures above 
~8000K, two times larger than the measured peak temperature. Auger heating (red arrows) 
dominates and leads to the formation of a single hot electronic distribution (b). In this case 
impact ionization and Auger heating are equally likely with zero net effect, and the electronic 
system can only cool down via phonon emission. 

conditions for carrier multiplication may be met for negligible doping of the graphene layer, 
smaller pump ffuences and higher excitation energies [8]. 

In summary, we have used time- and angle-resolved photoemission in the EUV as a 
tool to critically assess the potential of graphene for optoelectronic applications. Two key 
contributions are presented. 

Firstly, we show that population inversion occurs immediately after interband excita- 
tion, substantiating the view of graphene as a zero bandgap semiconductor. Previous time- 
resolved optical experiments had evidenced negative optical conductivity (optical gain) [3]. 
However, these measurements only probed at higher photon energies and could not detect 
the dynamics close to the Dirac point. The present experiment allows for a direct measure- 
ment down to arbitrarily small energies and with a temporal resolution that would not be 
possible for terahertz optical probes. Secondly, by directly measuring the band occupancy 
N(E,t) at all energies and times, we have made a quantitative evaluation of the relevance 
of carrier multiplication. 

Direct mapping of the non-equilibrium carrier distribution and primary accelera- 
tion/scattering events with few femtosecond or even attosecond resolution [35] is clearly 
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a highly interesting area of future research. Combining ultrafast photoemission tech- 
niques with nanoscale spatial resolution, for example in Photoemission Electron Microscopy 
(PEEM), would make it possible to investigate nanoscale carrier dynamics as well as real 
device transport. 

Methods 

Details about the sample growth are reported in Ref. [HJ [12]. The SiC(OOOl) substrate 
was etched in Hydrogen atmosphere to produce atomically flat terraces. In a second step, the 
substrate was annealed in inert Argon gas to grow a homogeneous carbon monolayer that 
was subsequently decoupled from the substrate by Hydrogen intercalation, resulting in a 
completely sp 2 -hybridized graphene layer. After growth, the samples have been transported 
in air and inserted into an ultra-high vacuum chamber, where the samples were annealed at 
about 500°C to remove adsorbates, and characterized with static ARPES. Measurements at 
the K-point perpendicular to the TK-direction were performed at the Max Planck Institute 
for Solid State Research in Stuttgart (Germany) with Helium II radiation at fko = 41 eV. The 
static ARPES measurement along the TK-direction was done at the Swiss Light Source (SLS) 
of the Paul Scherrer Institute (Villigen, Switzerland) with a photon energy of hu> = 30 eV 
and linearly polarized light. Measurements at the SLS were performed on a different sample 
[T2] (but grown after the same recipe) than the other static and time-resolved experiments 
in the manuscript. 

Time-resolved ARPES measurements were performed at the materials science end station 
at Artemis (Central Laser Facility, Didcot, United Kingdom), using synchronized near- 
infrared (NIR) and EUV photons as pump and probe pulses, respectively. The facility is 
equipped with a 1 kHz Ti: Sapphire laser system that delivers pulses with a central wavelength 
of 790 nm and a nominal pulse duration of 30 femtoseconds. NIR photons are generated 
via optical parametric amplification (OPA, 950 meV) and subsequent difference frequency 
generation (DFG, 300 me V). EUV photons in the range from 20 to 40 eV are obtained via 
high harmonics generation in an Argon gas jet. For the present experiment a photon energy 
of focjprobe — 31.5 eV has been chosen, using a time preserving monochromator to maintain 
a short pulse, enabling access to electrons emitted at the K-point of the two-dimensional 
graphene Brillouin zone. 
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All measurements in this paper were performed at room temperature, except the static 
dispersion along the TK-direction shown in Fig. 1(b) which was acquired at a sample 
temperature of 70 K. 
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